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Abstract 

The risk of extreme environmental events is of great importance for both the 
authorities and the insurance^] industry. This paper concerns risk measures in a 
spatial setting, in order to introduce the spatial features of damages stemming from 
environmental events into the measure of the risk. We develop a new concept 
of spatial risk measure, based on the spatially aggregated loss over the region of 
interest, and propose an adapted set of axioms for these spatial risk measures. 
These axioms quantify the sensitivity of the risk measure with respect to the space 
and are especially linked to spatial diversification. The proposed model for the cost 
underlying our definition of spatial risk measure involves applying a damage function 
to the environmental variable considered. We build and theoretically study concrete 
examples of spatial risk measures based on the indicator function of max-stable 
processes exceeding a given threshold. Some interpretations in terms of insurance 
are provided. 

Key words: Excursion sets; Extreme value theory; Max-stable processes; Mixing 
properties; Spatial dependence; Spatial diversification; Spatial risk measures. 


1 Introduction 


It is of prime importance for both authorities and insurance companies to take the spatial 
features of environmental risks into account. For authorities, it is crucial to be able to 
detect the areas at risk: is it safe to build houses in a given area or better somewhere else? 
Similarly, an insurance company has to choose its geographical zone of activity as well as 
its portfolio size. The last issue is obviously related to spatial diversification. Thus, tools 
(especially risk measures) capable of quantitatively dealing with spatial diversification 
are needed. 

The notion of risk measure has been widely studied in the literature. A risk measure 
II is a function from a set of random variables (typically a cone) to the real numbers, that 
satisfies some axioms. The seminal paper by Artzner et al. (1999) introduced the concept 


of coherent risk measure, which was then generalized to the convex case by IFollmer and 


Schied (2002) and Frittelli and Rosazza Gianin (2002). This static framework for risk 


measures was then extended to the conditional and the dynamic setting. For a detailed 
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review of conditional and dynamic risk measures, we refer to Acciaio and Penner (20 11 
The most widely used risk measure in the regulatory context is the Value-at-Risk (VaR) 
The aforementioned risk measures are univariate. In R, the natural order allows to 
easily define the notion of quantile and therefore VaR. However, in dimensions higher than 
1, the lack of such a natural order makes straightforward generalizations non-trivial. This 
is why many different definitions of multivariate quantiles have emerged in the literature. 


Di Bernardino (2013). 


For a detailed review of these, we refer to Serfling (2002). Regarding extensions of VaR 


to the multivariate setting, see e.g. Embrechts and Puccetti (2006) as well as Cousin and 


To the best of our knowledge, only Follmer (2014) and Follmer and Kluppelberg (2014) 


use the expression spatial risk measure. Consider that a local conditional risk assessment 
is carried out at each node of a network of financial institutions, in the sense that the 
risk measure applied takes into account the situation at the other nodes. The main issue 
they raise is whether the local risk assessments can be aggregated in a consistent way in 
order to provide a global risk measure. 

Since risk measures initially appeared to deal with financial risks, in an insurance 
context, they do not make explicit the influence of the region where the policies were 
underwritten. However, in an insurance portfolio, this particular region has an obvious 
impact on the risk undertaken by the company. In this paper, we introduce a new 
notion of spatial risk measure by explicitly disentangling the spatial region and the hazard 
generating losses over this region. Then, we study how the measure of the risk is expected 
to evolve with respect to some of the features of the spatial region, such as its location 
and its size. This leads to a set of axioms adapted to the spatial context. Contrary to 


the axioms proposed by Artzner et al. (1999), we study the sensitivity of the measure of 


the risk with respect to the space variable. 

Let us denote by A C R 2 the region under consideration and by {Cp(x)} xg if l 2 the 
process of the cost (in the whole paper, cost refers to an economic or insured cost) due 
to a particular environmental hazard (e.g. a hurricane). The easiest approach to build a 
spatial risk measure is to integrate a (univariate and static) risk measure (e.g. the variance 

1 f 

or VaR) over A, i.e. to consider —— / n(Cp(x)) dx, where |.| denotes the Lebesgue 

1^1 J A 

measure. However, if the process Cp is stationaryjq then the distribution of Cp{x ) is 
independent of x and thus the previous quantity is equal to n(C(0)). The corresponding 
spatial risk measure reduces to the univariate risk measure associated with a single site, 
meaning that this approach does not account for the spatial dependence structure of the 
cost process. In order to overcome this defect, we define our spatial risk measure by 
applying a (univariate and static) risk measure to the normalized aggregated loss over A. 

After having defined our general notion of spatial risk measure and the corresponding 
set of axioms, we introduce a general model for the cost process. The latter involves 
a mapping of the environmental variable under consideration to a cost via a damage 
function. Then, the paper focuses on the study of concrete examples of spatial risk 
measures based on a specific cost process involving the indicator function of threshold 
exceedances of the environmental variable. In a context of climate change, some extreme 
events tend to be more and more frequent; see e.g. Bevere and Mueller (2014). It is 


of prime importance for authorities as well as for the insurance industry to assess the 


2 Although it would be more rigorous to call this risk measure a quantile, we will mainly use the term 
VaR since it now entered in the common language of applications in banking and insurance, 
throughout the paper, stationarity refers to strict stationarity. 
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risk of natural disasters. A precise assessment of the risk of extreme events is crucial 
in order to satisfy capital requirements under the Solvency 11 regulatory framework. 
Therefore, due to the spatial feature of the environmental events, we model the process 
of the environmental variable using max-stable processes, which constitute an extension of 
extreme value theory to the level of stochastic processes (see e.g. de Haan, 1984; de Haan 


and Pickands, 

1986; 

Resnick, 

1987) 


The remainder of the paper is organized as follows. Our concept of spatial risk mea¬ 
sure, its corresponding set of axioms and a general model for the cost process are in¬ 
troduced in Section [2] Then, Section [3] describes the specific model for the cost process 
we use subsequently in order to build concrete examples of spatial risk measures. These 
examples are studied in Section |4j Section [5] concludes. Throughout the paper, the ele¬ 
ments belonging to R d for some d > 2 are denoted using bold symbols, whereas those in 
more general spaces will be denoted using normal font. All proofs can be found in the 
Appendix. 


2 Spatial risk measures 

Let A be the set of all compact subsets of R 2 with a positive Lebesgue measure. Denote 
by V a family of distributions of real-valued stochastic processes on R 2 having locally 
integrable sample paths. Each process represents the cost caused by the events belonging 
to specified classes and occurring during a given time period, say [0, Tfi . In the follow¬ 
ing, Tl is considered as fixed and does not appear anymore for the sake of notational 
parsimony. The events considered here have a spatial extent and thus it is natural to 
consider a cost process on R 2 . Each class of events (e.g. a heat wave or a hurricane) will 
be referred to as a hazard in the following. Let £n be the set of all real-valued random 
variables defined on an adequate probability space. A risk measure typically will be some 
function II : £n l— t R- 

2.1 Definitions 

We first give the definition of the normalized spatially aggregated loss, which allows to 
disentangle the contribution of the space and the contribution of the hazards. Indeed, in 
the case of an insurance company, the total loss in a portfolio of risks depends on both 
the region where the policies have been underwritten and the hazards covered in these 
policies. 

Definition 1 (Normalized spatially aggregated loss). For A 6 A and P G V , the nor¬ 
malized spatially aggregated loss is defined by 

Ln(A , P) = t^j C'p(x) dx, (1) 

where the stochastic process {Cp(x)} xg iR 2 has distribution P. 

The quantity L(A, P) = / C(p(x) dx is finite since A is compact and Cp has locally 

J A 

integrable sample paths. The random variable L(A, P) corresponds to the total economic 
or insured loss over region A due to specified hazards and is therefore of interest for 
spatial risk management. It seems more relevant, for both theoretical study and practical 
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interpretation, to consider the normalized spatially aggregated loss which is a loss per 
surface unit and can be interpreted in a discrete setting as the loss per insurance policy. 

Using the concept introduced in Definition [lj we now define our notion of spatial risk 
measure, which makes explicit the contribution of the space in the risk measurement. 

Definition 2 (Spatial risk measure). A spatial risk measure is a function 1Zn that assigns 
a real number to any region A E A and distribution P E V : 

Vn : A XV —y R 

(A,P) ^ n n (A,P) = U{L N (A,P)), 

where Ln(A,P) is defined in ([!]). 

If the distribution P of the cost process is given, then the function TZ\fi-,P) sum¬ 
marizes, for any region belonging to A, the risk caused by the hazards characterized by 
P. In the following, 7^n(-,P) will be referred to as the spatial risk measure induced 
by P. For many useful risk measures II (variance, VaR, Expected Shortfall, ...), this 
notion of spatial risk measure allows to take (at least) part of the spatial dependence 
structure of the process Cp into account. However, this is not true in the trivial case of 
the expectation due to its linearity. 

It now appears natural to analyze how Pn(Al, P ) evolves with respect to A for a given 
P. Some natural properties of 7Zn(-,P) are described in the set of axioms presented 
below. The spatial properties of (-,-P) depend on both the risk measure n and the 
probabilistic properties of the cost process characterized by P. 

2.2 A set of axioms for spatial risk measures 

This section provides a set of axioms in the context of the spatial risk measures introduced 
above. These axioms concern the spatial risk measures properties with respect to the 
space and not to the cost distribution, the latter being considered as given by the problem 
at hand. 

Definition 3 (Set of axioms for spatial risk measures). For A E A, let b a denote its 
barycenter. For a fixed P E V, we define the following axioms for the spatial risk measure 
induced by P: 

1. Spatial invariance under translation: 

for all v E R 2 and A E A, 1Zyi(A + v, P) = 7Zp(A, P), where A + v denotes the 
region A translated by the vector v. 

2. Spatial sub-additivity: 

for all Ai, A 2 E A, V^Ai U A 2 , P) < min{Pn(A? -P), 1Zyi(A2 1 P)}. 

3. Asymptotic spatial homogeneity of order —«, a > 0; 

for all convex A E A, 

K n (\A,P) + ^ + (2) 

where A A is the area obtained by applying to A a homothety with center and 
ratio A > 0, and K\ E R, K 2 E R\{0}. The constants A\ and K 2 can depend on 
A. 


4 


These axioms are well defined. For all A 1; A 2 compact, Ai U A 2 is compact as a finite 
union of compact sets. For all A e A, v e R 2 and A > 0, A + v and A A are compact as 
images of compact sets by continuous functions. 

ft is also possible to introduce the axiom of spatial anti-monotonicity: for all 
A \, A 2 e A, Ai C A 2 1Zu(A 2i P) < 1Zn(Ai, P). ft is easily shown that the latter is 
equivalent to the axiom of spatial sub-additivity. 

The axioms proposed in Definition [3] seem relevant when the cost process Cp is sta¬ 
tionary. In this case, there is clearly spatial invariance under translation. Our spatial 
sub-additivity axiom means that the risk associated with the normalized spatially aggre¬ 
gated loss is lower when considering the union of two regions instead of only one of these. 
It indicates that there is spatial diversification, which appears as a natural property when 
Cp is stationary. If this axiom is satisfied, an insurance company would be well advised 
to underwrite policies in both regions A\ and A 2 . As we will see in the following, the 
axiom of asymptotic spatial homogeneit}^] of order —a seems relevant, especially for II 
being the expectation, the variance or VaR. 

Although there are some links between our notion of spatial risk measures and financial 
risk measures as for instance summarized in Follmer and Schied (2004), the inclusion of 
the space and the C p process in Definition [2] sets our approach rather aside. 

For the study of concrete examples of spatial risk measures as well as practical appli¬ 
cations, we need to specify a model for the cost process Cp. A potential general model 
is provided in the following section, although the spatial risk measure framework defined 
above can be used in a much wider setting. 


2.3 A general model for the cost process 

The model we propose for the cost process {Cp(x )} xgR 2 requires two components. The 
first one concerns the cost generating hazards. We assume that the cost is only due to a 
unique class of events, i.e. to a unique natural hazard. The latter (e.g. a heat wave or 
a hurricane) is described by the stochastic process of an environmental variable (e.g. the 
temperature or the wind speed, respectively), denoted by {Z(x)} xgR 2 . We assume that 
Z is representative of the risk during the whole period [0, Tp\. 

The second component involves a model mapping the natural hazard into a damage 
and thus a cost. This model requires both the destruction percentage and the exposure at 
each location. The destruction percentage is obtained by applying a damage function (also 
referred to as vulnerability curve in the literature), denoted by D, to the natural hazard. 
This damage function is specific to the type of hazard considered. The exposure process, 
denoted by {A'(x)} xeR 2 , can be considered as deterministic and involves especially the 
demographic, economic and topographic conditions. Finally, the destruction percentage 
must be multiplied by the exposure, yielding the following model for the cost process: 


{Cp(x)} x ge = {E(x) £>(Z(x))} XER1 . 


(3) 


If insured losses are of interest, then the specific terms of the insurance policies should 
also be taken into account. 

Note that many generalizations of the model in (|3]) are possible, such as the intro¬ 
duction of a dependence of D with respect to the location in order to account for the 


4 Connected to the study of the function A H> 7\Ln(A A, P) is, in a temporal context, the analysis of the 


term structure of risk measures. For related studies, see e.g. Embrechts et al. (2005) and Guidolin and 


Timmermann 


(2006). 
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varying type of building in space, the explicit introduction of the temporal aspect and 
the consideration of several environmental variables (for instance, in case of hurricanes 
where damages depend on both the wind speed and the rainfall amount, or in case of 
multi-hazards). 

Remark 1. The presence of P in the right-hand term of (J3]) is implicit: the distribution 
P of the process Cp indeed depends on the three components of the right-hand term. 


3 A specific model for the cost process 

3.1 The model 

In order to study concrete examples of spatial risk measures in the following, we consider 
a particular case of the model for the cost introduced above. For the purpose of this 
paper, we choose the exposure to be uniformly equal to unity. Moreover, we consider, for 
a threshold u > 0, the damage function 


D{z) ^ > 0, 


which is adapted to the impact of high values of an environmental variable. In order to 
obtain a more realistic destruction percentage, we could for instance take 


D(z) 




k 


u 


1 


where k u > 1 depends on u. Nevertheless, there is no loss of generality in setting k u = 1. 

An important focus of our paper is the modeling of losses stemming from extreme 
events; this is particularly relevant for authorities as well as insurance companies. Hence, 
for all x G R 2 , we consider Z(x) to be the maximum of the considered environmental 
variable at location x over the period of interest [0, Tl\. This choice implies that the 


loss L(A,P) = / Cp(x) dx corresponds to the spatial aggregation of the costs caused 

J A 

by the worst events happening at each location during the time period considered. As 
explained in the next section, in this case, max-stable processes (see e.g. de Haan, 1984 


de Haan and Pickands 

1986, 

Resnick, 

1987) 


1987) are ideally suited for modeling purposes. 


Thus, for the remainder of the paper, we will assume the process {Z(x )} xgR 2 to be 
max-stable. Moreover, we will often make the classical assumption that Z has standard 
Frechet margins, i.e., for all x G R 2 and z > 0, 


P(Z(x) <z)= exp 



A max-stable process having standard Frechet margins will be referred to as a simple 
max-stable process. 

To summarize, the specific cost process we consider in the following is 


{Cp ( x )} xeR 2 - {l{z(x)> u }} xgR2 , (4) 

where {Z(x )} xgR 2 is a max-stable process and u > 0. Note that Cp has locally integrable 
sample paths since |C*p(x)| < 1, for all x G R 2 . 
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3.2 A short introduction to max-stable processes 

In the following, U \J V denotes the supremum and d G JM\{0}. Moreover, = and - 4 - stand 
for equality and convergence in distribution, respectively. 


3.2.1 Definition and comments 


Definition 4 (Max-stable process). A stochastic process {G(x)} xeRd is said to be max- 
stable if there exist sequences of continuous functions (ar(x),x G R d ) T>1 > 0 and 
(&t( x ), x ^ R d ) T> i G R such that if {Gt(x)} xgR d, t = 1,... ,T, are independent replica¬ 
tions of G, 


Vf = i(G f ( X )} - 6r(x) 
a T (x) 


— {G'(x)} xeR d. 


ceR d 


Now, let {Tj(x)} xgRt i, i — 1,... ,n, be independent replications of a stochastic process 
{T(x)} xgR d. Let (c n (x),x G R d ) n>1 > 0 and (d n (x),x G R d ) G R be sequences of 
continuous functions. If there exists a non-degenerate process {G'(x)} xgR d such that 


Vr=l(^( X )} - dn(x) 


xeR d 


->• {b'( x )}xeR d > for n —)• oo, 


(5) 


then G is necessarily max-stable, see e.g. de Haan (1984). Therefore, max-stable processes 


are well suited to model the joint behavior of the temporal maxima at all points in space. 

For practical purposes, the number n of observations over which the maxima are taken 
depends on the length T L of the time period considered. We assume that the limit in 
(|5| has been reached. This assumption of course needs some statistical justification in 
the examples analyzed. Classically, Tl — 1 year and n = 365. However, an insurance 
company can be more interested in the loss due to some particular events than in that 
corresponding to the worst event in the year. In that case, we can for instance take 
Tl — 1 week. 


3.2.2 Spectral representations and max-stable models 

Different spectral representations of max-stable processes can be found in the literature. 
Here, we present a very general one which can be found in Strokorb et al. (2015) and is 
especially based on de Haan (1984), Schlather (2002) and Kabluchko (2009). Let (H, 7F, v) 
be a probability space and let ipi denote the Lebesgue measure on the Borel c-algebra B d 
of R d . Moreover, let / : R d x hi t—)■ [0, oo] be a measurable function. 


Theorem 1. Simple max-stable processes that are separable in probability allow for a 
spectral representation of the form 


{Z(x)} x6Ei i | \/ {t/,l4(S,)} 1 , (6) 

v*=i J xeR d 

where the ( Ui , S))* are the points of a Poisson point process on (0, oo) x E with intensity 
u~ 2 du x p(ds) for some Polish measure space ( E , £, p) and the functions k x : L G (0, oo) 
are measurable such that f E V x (s) p(ds) = 1, for each x G R d . Moreover, any process of 
the form © is a simple max-stable process. 
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This general spectral representation especially includes two famous classes of repre¬ 
sentation: 


1. (Mixed) Moving Maxima representation 


Let / satisfy 



/(x,u;) dx h'(du) = 1, and let (E, £, p) = (R d xf2, v d xv). 


n Jn d 


Considering the following spectral functions in (J6]) , 

\4(s,ta) = /(x — s,cj), for all (s,cj) G H d x hi and x G R d , 

yields the so-called Mixed Moving Maxima processes (M3 processes) or Moving 
Maxima processes (M2 processes) if the random shape function is deterministic 
(i.e. if v charges only one point ujq G hi). 

2. Stochastic processes-based representation 

Let / satisfy / /(x, ui) v(duj) = 1, for all x G R d . We consider (E, £, p) = (hi, J-, v) 

Jn 

and V x (u) = /(x, cu), for all u G hi and x G R d . It is equivalent to denote 

14(.) = V(x), for all x G R d , 
where {V(x)} xeR d is a stochastic process. 

For specific cases included in Theorem [lj we refer the reader to de Haan and Ferreira 


(7) 


(2006), Corollary 9.4.5 and Theorem 9.6.7, Smith (1990), and Schlather (2002), Theorems 
1 and 2. 

The probabilistic representation (J6| has led to different parametric models for max- 
stable processes, some of which are presented in the following. Let {£(x)} xgE{ d be a 
stationary standard Gaussian process with any correlation function p. 

The Smith model 


Smith (1990) introduces a specific M2 process where the deterministic shape function / 


is the density of a d-variate Gaussian random vector with mean 0 and covariance matrix 
E. 

The Schlather model 

Schlather] (2002) proposes to set V(x) = a/27t e(x) in (0). All correlation functions 


stemming from the geostatistical literature can be used, allowing for a rich diversity of 
behaviors. 

The Brown-Resnick model 


Kabluchko et al. (2009) introduce a generalization of the model introduced by Brown and 


Resnick 


(1977), by taking in (0) V(x) = exp ( kF(x 




, where {fF(x)} xeR d is a 


Gaussian process with stationary increments and crj^(x) = Var(W(x)), for all x G R d , 
where Var stands for the variance. This process is referred to as the Brown-Resnick 
process and its distribution depends only on the semivariogram of W, defined by 7 (h) = 
|Var(W(x + h) — W(x)), h G R d . The special case where W(x) = cr £ e(x), with cr £ > 0, 
leads to the so-called geometric Gaussian process, see e.g. Davison et al. (2012). It should 
also be noted that the Smith process is a particular case of the Brown-Resnick process; 
see e.g. Yuen and Stoev (2013), Section 2.2. 

A correlation function p is said to be isotropic if, for all X!,x 2 G R d , p(x 1; x 2 ) only 
depends on ||xi — x 2 ||, where ||.|| denotes the Euclidean distance. In the cases of the 













































Family 

Correlation function (h > 0) 

Range of validity 

Whittle-Matern 

p( ft ) - I'm U) a «U) 

/ / \ 2\ -C2 

c 2 > 0 

Cauchy 

pW = (1 + ( A ) ) 

c 2 > 0 

Powered exponential 

P(h) = exp (- (A)") 

0 < c 2 < 2 


Table 1: Some families of isotropic correlation functions. Here, C\ and C 2 are the range 
and smoothing parameters, T is the Gamma function and K C2 is the modified Bessel 
function of the third kind of order C 2 . 


Schlather and geometric Gaussian processes, we will compare and contrast the families 
presented in Table [l] The powered exponential correlation function with C 2 = 1 will be 
referred to as the exponential correlation function. 


3.2.3 Extremal coefficient and mixing properties 

The extremal coefficient (see e.g. Schlather and Tawn, 2003) is a measure of spatial 
dependence for max-stable processes. Let {Z(x)} xgR d be a simple max-stable process. 
In the case of two locations, the extremal coefficient function 0 is defined by 


P 


(Z(x. 1 ) < u, Z (x 2 ) <u) = exp > x i , x 2 e R d 


If Z is stationary, then 0(x l5 x 2 ) only depends on h = x*—x 2 . Similarly to the correlation 
function, 0 is called isotropic if, for all xi,x 2 £ R d , 0(xi,x 2 ) only depends on ||xi — x 2 ||. 

Results about spatial diversihcation we present in Section [4] are expressed in terms of 
the extremal coefficient and it turns out that the latter is directly related to mixing prop¬ 
erties of max-stable processes. For all h £ R, let us define r(h) = 2 —0(h), where © is the 


extremal coefficient function. Kabluchko and Schlather (2010) prove (Theorems 3.1 and 


3.2, respectively) that a simple, stationary and measurable max-stable process {Z(x )} x£ ]a 

1 f l 

is mixing if and only if (iff) lim^oo r(h ) = 0 and is ergodic iff lim^oo - / r(h ) dh = 0. 

‘ Jh=l 


These results can be extended to R d ; see e.g. Dornbry (2012), p.20. 


Let us now consider the process {i/(x)} xgR d = {.D(Z(x))} xgR d, where D is a real¬ 
valued and measurable function defined (at least) on the state space of Z. It is worth 
mentioning that if Z is mixing, this also holds true for H. Indeed, the a-algebra generated 
by H is included in that generated by Z. This means that the mixing properties of the 
environmental process are also valid for the cost process Cp. 


4 Some concrete examples of spatial risk measures 

In this section, we consider the normalized spatially aggregated loss obtained by combin¬ 
ing (fi} and Q: 

1 f 

L N (A,P) = — jj.{z {x) >u} dx, (8) 

where {Z(x )} xgR 2 is a max-stable process and u > 0. Note that, for all A £ A and 
P £ V, | Ln{A, P)| < 1, meaning that all moments of Ln(A, P) are finite. This quantity is 
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interesting when analyzing, for instance, the impact of high temperatures on populations 
as well as on the distribution network of electricity, typically as in the case of the European 
heatwave in 2003. If u is well chosen, Ln{A,P) represents indeed the proportion of the 
surface area at which the temperature exceeds a dangerous threshold for populations and 
electric cables. 


It should be noted that / I{z( x )>m} dx corresponds to the so-called intrinsic volume of 

J A 

the excursion set E U (Z, A) = {x e A : Z(x) > u}. Excursion sets of stochastic processes 
have been widely studied in the literature; see e.g. Adler et al.| (2013), Spodarev (2014) 
and references therein. 

The dependence of L^(A, P) with respect to P lies in the distribution of Z. In the 
following, different max-stable models, especially those introduced in Section 3.2.2 will 


be considered. In each case, the model used will be explicitly indicated and therefore, we 
make the dependence in P implicit in Ln{A, P): from now on, Ln(A, P ) will be denoted 
Ln(A). Likewise, the spatial risk measure TZn(A,P) will be denoted IZn(A). 

The case of the expectation IZi(A) = E (Ljv(A)) is trivial if Z has identical margins. 
Using the linearity property of the expectation, we immediately show, if Z is simple, that, 
for all A E A, 


KM) 


1 — exp 



) 


meaning that 1Z\ does not depend on the region considered (and thus on its size). The 
spatial risk measure obviously satisfies the axioms presented in Definition [3j the order 
of the asymptotic spatial homogeneity being 0|^] The expectation is not a very useful risk 
measure since it does not involve any information relative to the variability. Furthermore, 
due to linearity, it does not account for the spatial dependence of the cost process Cp. In 
the following, we study the case of the variance in details, before providing some insights 
concerning VaR and briefly discussing some practical considerations. 


4.1 The variance 

We consider in this section the quantity IZ 2 (A) = Var (Ln{A)) . The variance allows 
taking into account part of the spatial dependence in the risk assessment. Hence, its 
study is interesting for both the risk management of extreme spatial events and the 
understanding of some properties of max-stable processes. As we will see, 7^2 is linked 
with the notions of spatial diversification, ergodicity and mixing. 

In the following, our aim is to study whether 1Z 2 satisfies the axioms presented in 
Definition [3] Before doing so, we study the function A 1 —> IZ 2 (AA) in detail. It is of 
course related to the property of asymptotic spatial homogeneity. However, we are also 
interested in the behavior of 7^2 (AA) for finite values of A. The latter can be of practical 
relevance for the insurance industry and, moreover, the results obtained will be used to 
prove the axiom of spatial sub-additivity. 

4.1.1 Study of IZ2 (AA), A > 0 

The expression of 1Z 2 (AA) in the case of simple max-stable processes is given in the next 
theorem. 

5 This is true for all cost processes having identical margins and finite expectation. 
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Theorem 2. Let {Z(x)} xgK2 be a simple max-stable process having 0 as extremal coef¬ 
ficient function. For all A G A and X > 0, we have 


iz 2 (xA) 


i 

AW 





— exp 



dx dy. 


(9) 


From Theorem [2j we can derive the expression of 1Z 2 {XA) in the case of simple max- 
stable processes having an isotropic extremal coefficient function, when A is either a disk 
or a square. Note that in the whole paper, disk and square refer to a closed disk with 
positive radius and a (closed) square with positive side, respectively. The result is given 
in the next corollary. 


Corollary 1. Let {Z(x )} xgIR 2 be a simple max-stable process having an isotropic extremal 
coefficient function 0(h), h > 0. Then: 


1. Let A be a disk with radius R. For all X > 0, we have 

r 2R 


IZfiXA) = — exp 


■ 2 i + 

u 


f d (h,R) exp 


'0 


0(A h) \ 
u 


where fd is the density of the Euclidean distance between two points 
and uniformly distributed on A, given by 


2h / 2 


f d (h, R) = — [ - arccos [ ^5 ) ~ —5\l 1 


7 r 


h 


2 R 


h 


7ri? 


h 2 

4K 2 


dh , (10) 

independently 


2. Let A be a square with side R. For all X > 0, we have 


1Z 2 (XA) = — exp 



f-Vzn 


+ 


f s (h,R) exp 




dh, 


( 11 ) 


where f s is the density of the Euclidean distance between two points independently 
and uniformly distributed on A, given by: 


• For h G [0, R], 

• For h G [R, R\f~2], 


fs(h,R) 


2vr h 8 h 2 2 h 3 

~W - ~B? + ~W' 


fs{h,R) 


^-2 - 6 + 3 V 6 - 1 + 


b+ 1 
Vfa1 


+ 2 arcsin 




2 h 

r r 


where b 


R 2 ' 


3. In both cases, if liiUA^-oo 0(A) exists, 
limiting risk measure given by 


then 1Z 2 (XA) converges as X ^ 00 to the 


— exp 



+ exp 


^ Hindoo 0(A) ^ 


( 12 ) 
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Remark 2. For more complex geometric shapes of the region A, little is known on the 
density of the distance between two points independently and uniformly distributed on A 


(see e.g. Moltchanov , 2012, Section 4.3.3), making this study more difficult. 


The following result directly follows from Corollary [lj 

Corollary 2. Let Z be as in Corollary [I] and A be a disk or a square. If, for all h > 
0 , 0 (h) = 1 (case of asymptotic perfect dependence whatever the distance), we have 

%(/!)= exp f-iCexpf-? 


Let us consider the max-stable models introduced in Section 3.2.2, in the case where 


they have an isotropic extremal coefficient function. This requires the correlation function 
p and the semivariogram 7 to be isotropic. Moreover, let us assume that p and 7 are 
non-increasing and non-decreasing, respectively. In this case, for h > 0, the function 
A 1 —>■ @(A h) is non-decreasing (see below), giving that A t— y IZ 2 (AA) is non-increasing for 
A being a disk or a square. Consequently, there is spatial diversification and H 2 (AA) 
has a limit as A —> 00 . Corollary [l] may be of interest for the insurance industry since 
it allows determining the characteristic dimension of the geographical area required to 
reach a specified low variance level. As shown in the following corollary, in some cases, 
the diversification can be total, in the sense that lim^oo TZ 2 (XA) = 0. 

Corollary 3. Let Z be as in Corollary^and A be a disk or a square. If lim,\->oo 0(A) = 2 
(case of asymptotic independence when the distance tends towards infinity), we have 


lim 72-2(A A) = 0 . 

A—^00 

Moreover, if lim^oo 0(A) exists, then liniA->.oc ©(A) = 2 iff lim^oo IZ 2 {\A) = 0. 


(13) 


This result is not surprising. Indeed, as mentioned in Section 3.2.3[ if additionally 
Z is measurable and stationary, liniA-^oo ©(A) = 2 implies that Z is mixing. Thus, 


the process I {z>uj is also mixing and therefore mean-ergodic, which implies (13) for all 
A E A. Furthermore, from Corollary [3j we deduce, as can be done from Kabluchko and 
Schlather (2010) for max-stable processes on R, that in the case of a simple, stationary 


and measurable max-stable process having an isotropic extremal coefficient function, if 
Hhia-hx) ©(A) exists, then ergodicity and mixing are equivalent. Indeed, if Z is ergodic, 
then I {z>u} is also ergodic and thus mean-ergodic, which implies that limA_).oo 7£ 2 (AA) = 
0 . Thus, if liniA-s.00 ©(A) exists, this gives that limA^.00 ©(A) = 2 , i.e. that Z is mixing. 

Theorem |2] shows that the behavior of the function A t —y 1Z 2 (\A) is mainly driven 
by the extremal coefficient function 0 ; the latter naturally depends on the max-stable 
model under consideration. In the following, we study the influence of A for some classical 
simple max-stable processes having an isotropic extremal coefficient function. Since the 
integrals in Corollary [T| have no closed form, we use a Riemann approximation. Without 
loss of generality, we set R = 1. The choice of the threshold u has no influence on the 
shape of the function A e->- IZfiXA) and we choose u — 1. 

The Smith model 


The extremal coefficient function is given by (see e.g. Davison et al. 2012) 


0(xi,x 2 ) = 2$ 


a/(xi — X 2 )'S X (X! X 2 ) 


xi,x 2 e R 2 , 
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where $ denotes the distribution function of a standard Gaussian random variable and 
the prime means transposition. The function 0 is isotropic iff £ is proportional to the 
identity matrix. Without loss of generality, let £ be the identity matrix. In this case, 
we have 0(h) = 2$ (|), for all h > 0. Hence, limA-^oo 0(A) = 2 and Corollary [ 3 ] gives 
that lim^oo IZ 2 (AH) = 0, meaning that the spatial diversification is total. We observe in 
Figure [I] that ^(AA) rapidly decreases to the limiting risk measure when A increases. 



Figure 1: The solid (respectively dotted) line depicts the evolution of ^(AH) with respect 
to A in the case of the Smith model, where A is a disk (respectively a square). The dashed 
line represents the limiting risk measure. 


The Schlather model 


The extremal coefficient function is given by (see e.g. Davison et al. 2012) 


1 - p(xi - x 2 ) 


0(x 1 ,x 2 ) = l + \l -— Xi ,x 2 gR 2 . 


The function 0 is isotropic iff the correlation function p is isotropic. In thi s case, 0 
1+ \J —1^-, for all h > 0. Moreover, for all h > 0, p{h) > — | (see e.g 


Abrahamsen, 1997 

Theorem 3.10), implying that, for all A > 0, 0(A) < 1.87, and thus limA-^oo77.2(XA 


(h) = 


> 0 . 


This shows that the process I {z> u } is not mean-ergodic, which is consistent with the 
fact that the Schlather process is neither mixing nor ergodic. In terms of insurance, this 
result means that the spatial diversification is never total. In the case of the correlation 
functions in Table[ll we have limA^-oo p{ A) = 0, giving, using (12), that limA-;.oo 77.2(AA) = 

— exp (—|) + exp ' 1+A 


We set the range parameter c\ = 1 and the smoothing parameter c 2 = 0.5. Figure [2] 
shows that the speed of decrease of 77 2 (AA) to the limiting risk measure depends on the 
type of correlation function. This decrease is much slower in the cases of the Cauchy and 
the powered exponential functions. The smoothing parameter also plays an important 
role, as can be seen from the comparison between the powered exponential case and 
the exponential one. Hence, the Schlather process allows for a large variety of spatial 
diversification behavior. 
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Figure 2: Each panel corresponds to a different correlation function. The solid (respec¬ 
tively dotted) line depicts the evolution of TZ 2 (A A) with respect to A in the case of the 
Schlather model, where A is a disk (respectively a square). The dashed line represents 
the limiting risk measure. 


On the whole, we observe that the decrease to the limiting risk measure is slower 
than in the case of the Smith model. This comparison must be made with equivalent 
characteristic distances of spatial dependence, meaning that the eigenvalues of £ must be 
close to Ci, which is the case. Note that the decrease is obviously slower when increasing 
the range parameter C\. 


The geometric Gaussian model 


The extremal coefficient function is given by (see e.g. Davison et al. 2012) 


0 (xi,x 2 ) 


2 $ 


'<?%[ 1 -p(x i x 2 ) 


xi,x 2 G R 2 . 


The function O is isotropic iff p is isotropic. In this case, 0(h) = 2<f> 
for all h > 0. Therefore, as before, lim^oo TZ^XA) > 0. In the cases considered here, 


a^[l-p(h)\ 
2 


liniA-ioc p( A) = 0, so that by ([12]), lim^oo R 2 (XA) = - exp (-|) + exp 


2*1 \ r - f - 


We set <7 e = 1 and, as previously, C\ — 1 and c 2 = 0.5. From the plots (available upon 
request), we can draw very similar conclusions to those related to the Schlather model. 
The main difference consists in the value of the limiting risk measure. 

In order to allow for faster spatial diversification, we consider a new instance of a M2 
process, the so-called tube model, defined below. 


The tube model 


Definition 5 (Tube process). The tube process is a M2 process where the deterministic 
shape function is written /( y) = hb I{||y||<J2 b } , y £ R 2 , with Rb > 0 and hb = -hj. 

o 

The density / is a tube with height hb and radius Rb centered at point 0. This process 
is simple and its bivariate distribution function is given in the next proposition. 
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Proposition 1. Leix 1 ,x 2 G R 2 , h = ||xi — x 2 1| andz\,z 2 > 0. The bivariate distribution 
function of the tube process is given by 


-log(P(Z(xi) < Z]_, Z (x 2 ) < z 2 )) 


where 


^[rtRl - A int (h )] + A if z 2 <z ll 
\ + ^[kRI - A int {h )] if Z2>z 1} 


A-int (^) 


2 (rI arcsin ~ ) “ f ~ */ h < 2i?b, 

0 if h> 2R b . 


The extremal coefficient function, which is isotropic, is given in the next corollary and 
depicted in Figure [3] for R b — 1. 


Corollary 4. The extremal coefficient function of the tube process is given by 


0 (h) 


1 - h h 


Rl arcsin 


y/AR 2 b - h 2 \ h 


2 Rb 




if 0 <h< 2 R b , 
if h > 2 R b . 



h 


Figure 3: Extremal coefficient function of the tube process. 


We can easily show that for all h > 0, the function A t— y 6(Xh) is increasing for A < 


2 R b 
h 


and constant above that. Thus, (10) and (JTTj) yield that the function A t— > TZ 2 (^A) is 
decreasing for A being a disk or a square. Hence, there is spatial diversification. An 
interesting property stems from the fact that 0 (h) reaches 2 whenever h > 2 R b , meaning 
that there is spatial independence at a finite distance. This explains intuitively why the 
spatial diversification is faster than in the case of the previously introduced models, as 
can be observed in Figure [5] (obtained with R b = 1). Obviously, lim^oo 0(A) = 2, and 
Corollary [3] yields linr^oo 1Z 2 (\A) = 0. Moreover, due to the spatial independence at 
finite distance, we have the following proposition. 
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Figure 4: The solid (respectively dotted) line depicts the evolution of 7^ 2 (AAl) with respect 
to A in the case of the tube process, where A is a disk (respectively a square). The dashed 
line represents the limiting risk measure. 

Proposition 2. Let {Z(x)} xg R 2 be the tube process and A e A. We have, for all A > 1, 
liniR 6 _).o 7^-2 (AA) = TZ 2 (A) = 0. 

The limit process arising as Rb tends to 0 corresponds to the case of asymptotic indepen¬ 
dence whatever the distance. 

By way of summary, a comparison of the four models considered (Smith, Schlather, 
geometric Gaussian and tube) is provided in Figure [5] In the cases of the Schlather and the 
geometric Gaussian models, a Cauchy correlation function has been used. These different 


CO 

co Smith Schlather Geometric Gaussian Tube 



Figure 5: Each panel corresponds to a different max-stable model. The solid (respec¬ 
tively dotted) line depicts the evolution of 7\1 2 (AA) with respect to A, where A is a disk 
(respectively a square). The dashed line represents the limiting risk measure. 

processes show a large variety of behavior, both in terms of speed and “completeness” of 
spatial diversification. The Brown-Resnick model itself includes various types of behavior. 
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The two particular cases considered here (the Smith model and the geometric Gaussian 
model) are very different. In all cases, the behavior of the function A i —> IZ 2 (AA) is very 
similar in the case of the disk and the square. The only difference consists in the fact that 
the spatial diversification is slightly slower for the square. For an insurance company, the 
characteristic dimension of the geographical area required to reach a specified low variance 
level depends to a large extent on the type of max-stable process driving extreme events 
and on the geometric shape of A. 

Remark 3. It is important to note that our analysis has been carried out with a stan¬ 
dardized characteristic distance of spatial dependence: the eigenvalues of E, the range 
parameter c\ and the radius of the tube Rb are equal to 1. However, in a specific case 
study, the characteristic dimension of the area required to reach a given level of variance 
depends on the real values of E, c\ or Rb on the region of interest. 


4.1.2 Axioms 

The following theorem shows that the spatial risk measure IZ 2 satisfies, at least partially, 
the different axioms introduced in Definition [3] under some conditions. 


Theorem 3. 1. For any stationary max-stable process {Z(x)} xg R 2 , the spatial risk mea¬ 
sure IZ 2 satisfies the axiom of spatial invariance under translation. In particular, this is 
true for the Smith, Schlather, Brown-Resnick, geometric Gaussian and tube processes. 

2. For any simple and stationary max-stable process {Z(x )} xgE 2 having an isotropic 
and non-decreasing extremal coefficient function, IZ 2 satisfies the axiom of spatial sub¬ 
additivity when the two regions are both a disk or a square. In particular, this is true for 
the Smith process having a covariance matrix proportional to the identity, the Schlather 
and geometric Gaussian processes having an isotropic and non-increasing correlation 
function, the Brown-Resnick process having an isotropic and non-decreasing semivari- 
ogram and the tube process. 

3. Let {Z(x)} xg R 2 be a simple max-stable process having an extremal coefficient function 
0 such that for all x, y £ R 2 , @(x, y) only depends on x — y. Thus, 0(x, y) is written 
B(x-y). If 

0 < / (2 — 0(x)) dx < 00 , (14) 

Jr 2 


then the spatial risk measure IZ 2 satisfies the axiom of asymptotic spatial homogeneity of 
order —2, with 


K\ = 0 and K 2 


a 


2 


\A\’ 


where, for u > 0, 


a = 



(15) 


In particular, this is true for the Smith process, the Brown-Resnick process with a semi- 
variogram satisfying 7 (h) = 77 1 1 h 11 “ 7 h £ R 2 , where rj > 0 and a £ (0,2], and the tube 
process. 

Remark 4. It turns out that the axiom of asymptotic spatial homogeneity of order 
—2 of spatial risk measures associated with the variance is satisfied for more general 
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cost processes than the one considered in this section. Indeed, let Cov denote the co- 
variance and {C'p(x)} xg R 2 be a cost process having locally integrable sample paths and 
such that, for all x,y G R 2 , Cov(Cp(x),Cp( y)) only depends on x — y, and a 2 g = 


Cov(Cp(0),Cp(x.)) fix G (0, oo). From the proof of Theorem 


R 2 


3. Bullet 3, we eas¬ 


ily see that the quantity Var 



CpCx) fix 


homogeneity of order —2 with K\ = 0 and It 2 = 


satisfies the axiom of asymptotic spatial 



4.2 Central limit theorem 


Having shown that, under some conditions, H 2 satisfies the axiom of asymptotic spatial 
homogeneity of order —2, the question whether Ln(XA) (for A being a convex in *4.) 
satisfies a central limit theorem arises pretty naturally. Using a result by Spodarev 
(2014), we can derive the following result. 


Theorem 4. Let {Z(x)} xg R 2 be a simple, stationary and measurable max-stable process 
having 0 as extremal coefficient function. If 



(2 — @(x)) fix < 00 , 


(16) 


we have, for all convex A E A, that 

A (L n (XA) -m) 4 - N ^0, ^ , for A -)• 00 , (17) 

where 

m = 1 — exp 

In particular, this is true for the Smith process, the Brown-Resnick process with a semi- 
variogram satisfying 7(h) = 7||h|| a , h G R 2 , where g > 0 and a G (0,2], and the tube 
process. 



It should be noted that there may be mixing max-stable processes for which (16) is 
not satisfied. This is the case if @(x) does not converge fast enough to 2 as ||x|| tends to 
infinity. 


4.3 The Value-at-Risk 


Here we focus on 77 3)Q (H) = VaR a (L N (A)), a G (0,1), where, for a random variable X 
with distribution function F, VaR„(A") = infja: G R : F(x) > a}. The risk measure VaR 
is of interest for insurance companies, and this mainly from a regulatory point of view 
(see e.g. McNeil et al. 2015). We adopt an approach similar to the case of the variance, 
giving first some insights about the function A 1 —> IZ^ )a {XA) before discussing whether 
lZz )0l satisfies the axioms of Definition [3} 
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4.3.1 Study of TZ 3t0l (XA), A > 0 

Since deriving a tractable formula for the VaR of L^{\A) is very difficult, in the case 
of max-stable processes which can be simulated, we propose to evaluate it as follows. 
We simulate the process under consideration on a grid containing the locations x m G 
A A,m — 1,..., A 2 M, where M is the number of sites in region A. Then, one realization 
of Ln(XA) can be approximated by a Riemann sum. By generating a number S of 
independent approximated replications of the random variable Ln(XA), an approximation 
of its distribution function is obtained. Finally, taking the empirical quantile of this 
distribution at level a yields an approximation of 1Z 3 , a - 

To illustrate this procedure, the evolution of 72-3,0.9(A A) with respect to A in the case 
of the Smith model is depicted in Figure [ 6 j for M — 49 and S = 10000. The sites are 
located on a regular grid. Of course, the larger M and S, the smoother the curve. Figure 



Figure 6 : 72-3,0.9(AR) with respect to A in the case of the Smith process, where A is a 
square of side 1 . 

7 seems to indicate spatial diversification. 


4.3.2 Axioms 

Although we do not have an explicit formula for 1Z 3 oi (\A), by taking advantage of The¬ 
orem [dj we know the asymptotic behavior (when A —y oo) of lZ 3}a (\A) under some 
conditions. This result is part of the following theorem. 


Theorem 5. 1. For any stationary max-stable process {Z(x)} xgR 2 ; the spatial risk mea¬ 
sure 723 jQ satisfies the axiom of spatial invariance under translation. In particular, this 
is true for the Smith, Schlather, Brown-Re snick, geometric Gaussian and tube processes. 
2. For any simple, stationary and measurable max-stable proces {Z(x )} xgR 2 with extremal 
coefficient function 0 such that (16) is satisfied, the spatial risk measure TZ 3 ^ a satisfies, 
for a 7 ^ 0.5, the axiom of asymptotic spatial homogeneity of order -1, with 


K\ = 1 — exp — 
u 


and Ko = 


O' gq 

•M' 


19 





where q a is the quantile at the level a of the standard Gaussian distribution. In partic¬ 
ular, this is true for the Smith process, the Brown-Resnick process with a semivariogram 
satisfying 7(h) = ?7||h|| a ? h € R 2 , where p > 0 and a £ ( 0 , 2], and the tube process. 

Since we do not have an explicit formula for H 3t0l (XA) for A < 00 , there is nothing 
we can say about spatial sub-additivity. However, the previous numerical study suggests 
that, in the case of the Smith model, the function A t —y lZ 3)a (XA) is non-increasing for A 
being a square, which would imply spatial sub-additivity when the two regions are both 
a square. 

Remark 5. From the proof of Theorem [5), Bullet 2, we see that asymptotic spatial ho¬ 
mogeneity of order — 1 of spatial risk measures associated with VaR is satisfied as soon 
as there is a central limit theorem for the normalized spatially aggregated loss. The exis¬ 
tence of such a theorem can be linked to the mixing properties of the cost process (see e.g. 
Theorem 1). 
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4.4 Practical considerations 


In our analysis, we have most often considered max-stable processes with standard Frechet 
margins although in practice, the univariate marginal distributions of max-stable pro¬ 
cesses are Generalized Extreme Value distributions, involving three parameters: p £ R, 
cr > 0 and £ £ R. These distributions are denoted GEV(/q a, £). For a random variable 

Y, we have that Y ~ GEV(/i, a, £) iff Z = [l + £ ? ~ standard Frechet. Hence, it 

is easy to show that, for u\ £ R, Y > u 3 iff Z > u, where 


u = 


[!+?(“)]' 

exp (gy) 


if £ 7^ 0, 

if £ = 0. 


(18) 


The parameters p, a and £ are such that u > 0. This means that the general case of 
identical GEV margins reduces to the standard Frechet case we have studied. A concrete 
application of the tools proposed in this section involves, in the case of identical margins, 
the following steps: 


1. Specify the threshold u\ of interest, for instance U\ = 30°G. 

2. Jointly fit to historical data the marginal parameters p, a and £ as well as the 
dependence parameters for different models (e.g. Smith, Schlather, Brown-Resnick, 
geometric Gaussian, tube), using, for instance, the composite likelihood approach 
(see e.g. 


Padoan et al. 2010 


3. Carry out model selection based, for instance, on the Akaike Information Criterion 
(Akaike 1974) or the likelihood ratio statistic (Davison, 2003). 


4. 

5. 


Transform u 1 in u by plugging the estimates of p, a and £ in (18). 


Study different spatial risk measures as explained, considering the threshold u com¬ 
puted in the previous step. 


In some situations, the marginal parameters may vary across space, meaning, using (18), 
that the threshold also varies across space and is now of the form w(x), x £ R 2 . In this 
case, 77-2(A A) involves the full expression of the bivariate distribution function and not 
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only the extremal coefficient, leading to more complex (but still numerically computable) 
expressions than <§• ( [To] ) and ( [n| ) . If the dependence parameters also vary across space, 
the complexity becomes even higher. 


5 Conclusion 

The first part of the paper introduces a new notion of spatial risk measure, based on the 
normalized spatially aggregated loss, and proposes a set of axioms adapted to the spatial 
context. These axioms appear natural for stationary cost processes. Contrary to the 
classical literature relating to risk measures, our axiomatic approach aims at quantifying 
the sensitivity of the risk measurement with respect to the space variable. Our motivation 
comes from insurance applications where losses are due to hazards having a spatial extent, 
typically such as environmental hazards. In order to develop concrete examples of spatial 
risk measures, we propose a general model that maps the process of the environmental 
variable into a cost, via a damage function. 

The second part of the paper considers a specific model for the cost, where the en¬ 
vironmental process is max-stable and the damage function is the indicator of threshold 
exceedances. Some spatial risk measures based on this model are studied from a theo¬ 
retical viewpoint. Under some conditions mainly on the underlying max-stable process, 
they satisfy, at least partially, the axioms proposed. This feature is true for more general 
cost processes, as shown especially in Remarks [4] and [5} Characterizing all univariate risk 
measures and processes allowing the proposed axioms to be satisfied would be relevant 
and useful. 

In the case of non-stationary cost processes, more adapted axioms should be proposed, 
depending of course on the type of non-stationarity. If the cost process can be assumed 
to be piecewise stationary, the concepts introduced here can be independently applied to 
each sub-region on which stationarity may be assumed. Asymptotic spatial homogeneity 
naturally requires sufficiently large sub-regions. 
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A Appendix: Proofs 

A.l For Theorem [2] 

Proof. We have, for all A G A, that 


E ([L(A)] 2 ) = E f (j I {z(x|> „, dx) J = E ( 


= E 



I{z(x)> u } dx J I{z{ y )>u\ dy'j 

z(~x.)>u,z(y)>u } dx dy 


A J A 



P(Z(x) > u, Z(y) > u ) dx dy. 


A J A 


Moreover, using the fact that Z has standard Frechet margins, we have 

P(Z(x) > u,Z(y) > u) — 1 + P(Z(x) < u, Z{ y) < u) — 2 exp (— — j , yielding 


E ((i(^4)) 2 ) = \ Af 


Hence, 


n 2 (A) = 


H 2 

-|H | 2 


lA 


1 - 2 exp — 
u 


1 - 2 exp — 
u 



+ I j ex p l v ~ <5> ^’ y ^ > ) d y- 


+ I j ex p l v ~ <5> ^’ y ^ > ) dx d y 



IA 2 



1 — exp — 
u 


. ©( x >yA ( 2 

exp-— exp — 

u 1 \ u 


A J A 


dx dy. 


The result is obtained by replacing A with AH. 


□ 


A.2 For Corollary [T] 

Proof. 1. Let A be a disk with radius R. For any integrable function g : (x, y) G 
R 2 x R 2 i-> R depending only on h = ||x — y||, it is easy to show that, for all A > 0, 


r> 2 A R 


/ / gfx,y) dx dy = A 4 |H|- / f d (h,XR) g(h ) dh. 

'\A JXA Jo 


(19) 


Moreover, we have (see e.g. Moltchanov 2012, Section 4.3.1) that 


fd{h , R) = | Fa.ccos (A) -f-Jl-fd | . for 0 < . < 


yielding, for all A > 0, f d (Xh,XR) = — f d (h,R). Therefore, by the change of variable 

h . 

rid = y, we obtain 
A 


p 2 XR /* 2 R 

/ f d (h , A R) g(h ) dh= f d (h d , R) g(Xh d ) d/i d . 

'o do 


( 20 ) 
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Applying (19) and (20) to g(x, y) = exp 


e(x, y ) 


u 


is 


assumed to be isotropic), and using , we obtain 


= exp 


e(||x-y| 


u 


(since 0 


TZ 2 {\A) = - exp /d(^> R ) ex P 


(-«) 


dh. 


2. The same reasoning yields the result for A being a square. The distribution function 
of the Euclidean distance between two points independently and uniformly distributed 
on a square with side R is written (see e.g. Moltchanov 2012, Section 4.3.2) 


F s (h, R) 

( n h 2 8 h 3 h A 

I If ~ FF n + w 4 

\ i b 2 2 

- — 2b — — + - a Jib — l) 3 + 2 \fb — 1 + 2b\/b — 1 + 2b arcsin 

V o Z o 



if h G [0, R], 
if h e [R,RV 2 ], 


where b 


F 

R 2 ' 


Hence, for h G [0, R], the density is 


fs{h,R) 


2tt}i 8 h 2 2b? 

If ~ W + If' 


For h G [A, Ry/2\ , we obtain 


fs{h,R) 


^-2 - b + 3y/b - 1 + 


b + 1 


+ 2 arcsin 


2-b 



2 h 
R 2 ' 


3. Since fd is continuous on the compact set [0, 2 R\, it is bounded. Thus, using Lebesgue’s 
dominated convergence Theorem and the fact that fd is a density, we obtain 


lim [ fd(h,R ) exp (— —— dh 

A->°o J g \ U J 


c2R 


exp 


fd(h,R ) exp 
lim A ^ 0O 0(A) 


lim 


A—^oo 


U 


e(A) \ 


u 


dh 


which gives the result when A is a disk. The argument is the same in the case of the 
square. □ 
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A.3 For Proposition [T] 

Proof. Since the (f/j,SA are the points of a Poisson point process on (0,oo) x R 2 with 
intensity u~ 2 du x ds, we have, for all z i, z 2 > 0 , that 


- log(P(Z(xi) < z l7 Z(x 2 ) < z 2 )) 


I R 2 


Zl 


Z2 


u z du ds 


/(xi - s)’/(x 2 - s 
f /(xi - s) /(x 2 - s) 

R 2 l *1 -2 

f /(x 1 - s' 


= / max 


ds 


-I 


'R 2 ^1 

r M h 


-a) /( x 2 ~ s ) ] dS + / 
z 2 / Jr 2 


f /(Xj-s) /(x 2 -s) 1 

l 2 1 2 2 J 


/(x 2 


^2 


If /(xj-s) > /(x 1 -s) 1 

l 2 2 ~ 2 1 / 


ds 


/(a) v. /( s + x 2 —■ 




ds -j- 


/• /(s 

/ - J 1 f/M> 

'R2 2^2 l *2 “ 


If f (s! /(s+x 1 —X 2 ) 1 


1 


J R 2 l 2 l ' 2 2 J JR 2 *2 l 2 2 - 2 l 

I P fill > /(X + ^-xQ A + j_ p / /(X) > /(X + x, - x 2 ) 

2 1 \ 2 1 2 2 / 2 2 V 2 2 2 1 


( 21 ) 


where X is a random vector having density /. 

r/(x). /(x+x 2 — Xl )i r _ 

“W - > -~-| - i/2 I{||X||<R 6 } > -1 I{||X+x 2 -xi||<R 6 })- 


Let Ei be the event 


V 


If z 2 > z u Ei = {||X|| < Rb}, giving P(H] ) = 1. Indeed X has density / and then 
||X|| < R b almost surely (a.s.). If Z\ > z 2 , E\ = {||X|| < R b and ||X + x 2 — xi|| > R b } = 
{||X + x 2 — xx|| > R b } since ||X|| < R b is satisfied a.s. Therefore, 


P(£i) 



I{||x+x 2 —xi||>f4} /( x ) dx 


hb / I{||x||<* 6 n ||x-(xi-x 2 )||>f?b} dx 

J R 2 

h b [tcRI - A int (h)\, 


where A int {h ) is the area of the intersection between the base of the tube of center 0 and 
that of the tube of center (xx — x 2 ). Note that by symmetry, the area of the intersection 
between the base of the tube of center 0 and that of the tube of center (x 2 — xi) is also 
equal to A int (h). 


J /(X) /(X + xi-xo)! _ 


- {Zl I{||X||<f4} > ^2 I{|jX+xi—x 2 ||<Ri,} } • 


Let E 2 be the event -[ ——- > 

l ^2 Zi 

If z 2 > zi, E 2 = {||X + xi - x 2 1| > R b }, giving P(-E 2 ) = h b [nRl - A int (h)\. If zi > z 2 , 
P (E 2 ) = 1 since ||X|| < R b a.s. 

Hence, using (21), we obtain 


- log(P(Z(xi) < Zi, Z(x 2 ) < z 2 )) 


%[nRb~ A int (h)} + A if z 2 < zi, 
A + ^[ 7 tRI - A int (h)\ if z 2 > zi. 


Since we consider the Euclidean norm, the bases of the tubes are circular. Let us 
then compute the area of the intersection between two discs, respectively with radius R b 
and centers Ci and C 2 at a distance h. This intersection is not empty iff h < 2 R b . We 
consider this case, represented in the following picture: 
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Let us introduce p = |(2 R b + h). By using Heron’s formula, the area of the triangle 
IC 1 C 2 , denoted by At, is given by 


A t = y/p(p - Rb)(p - Rb)(p - h). 


(23) 


hH 


Furthermore, denoting H the height of the triangle IC 1 C 2 , we have At = , giving 


H = 


2Ai 


h 


(24) 


H 


Denote a and f3 the angles IC 1 C 2 and IC 2 C 1 , respectively. We have sin a = sin (3 = —, 

Rb 


yielding, using (24), 


a = (3 = arcsin 


2 A t 

hRb 


(25) 


Let S be the area of the angular sectors delimited respectively by the angles a and f3. 
We have 

aR l (26) 


s = 


Combining ( |23[ ) , ( |25p and ( J26]) , we obtain, for h < 2 R b , 

2 Vp(p~ R bY(P~ h ) 


Aint(h) = 2(2 S — A t ) = 2 \ Rl arcsin 


hRb 


- y/p[p - R b ) 2 (p - h ) 


We finally obtain 


Aintijl) 


2 Rl arcsin 


( y/iRl -h* 




2 R b 


) - - hA if A 


< 2 / 6 ,. 

if h > 2 R h . 


The combination of (22) and (27) yields the result. 


(27) 


□ 


A.4 For Corollary [4] 

Proof. By definition of the extremal coefficient function, for all Xi,X 2 e R 2 , 0(xi,X2) = 


— log(P(Z(xi) < u, Z(x 2 ) < u) u. Hence, using (22), we obtain 


0(x!,x 2 ) = u [irR 2 b - A int {h)\ + ^ =2 - h b A int (h). 


The combination of (27) and (28) gives the result. 


(28) 

.O 
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A.5 For Proposition [2] 

Proof. We consider the case of A being a disk; the proof is exactly the same in the case 


of a square. Using (10), we have 


72-2 (A A) = - exp (~f^j + j fd(h,R ) exp 

t-2 R 


©(Ah) A 
u 


dh 


'2 R b 


fd(h,R ) exp ( dh. 


u 


Moreover, for A > 1 and h > 2Rb, we have @(A h) = 0(h) = 2. Hence, for all A > 1, 

ft 2 (A A) 

= [ ^ f d (h,R) exp dh — exp (--) + [ ^ f d (h,R) exp (dh 


u 

f d (h,R) exp dh 

r2Rb r / e(Ah)\ 


u 


u J 


= n 2 (A)+ / f d {h,R) 
Jo 


exp 


u J 


exp 


e(h) A 

u ) 


dh. 


When Rb tends to 0, the second term vanishes, giving lim^^o 7^ 2 (AH) = 1Z 2 (^4) • More¬ 
over, if Rb = 0, for all h > 0, for all A > 0, 0(Ah) = 2. Thus, (10) gives that for all 
A > 0,7^2(AA) = 0. □ 


A.6 For Theorem [3] 

Proof. 1. If Z is stationary, it is also true for the process Cp = I {z>u}- All processes 
mentioned are stationary. See e.g. Schlather (2002), Theorem 2, for the Schlather process 


and Kabluchko et al. (2009), Theorem 2, for the Brown-Resnick process; the tube process 


is obviously stationary as an instance of M2 process. 

2. Let us consider two regions A\ and A 2 being both a disk or a square and such that 
A\ C A 2 . Due to spatial invariance under translation, the region A 2 can be translated 
to region A' 2l where A' 2 corresponds to the region obtained by an homothety of Ai, with 
center the center of Ai and ratio denoted A > 1. Thus, TZ 2 (A 2 ) = TZ 2 (A' 2 ) = T^AA)- 
Moreover, for all h > 0, the function A i—> 9(Xh) is non-decreasing. Thus, (10) and 0 
give that 1Z 2 (\A) is a non-increasing function of A. Hence, we have 1Z 2 (\Ai) < 1Z 2 (Ai), 
giving 1Z 2 (A 2 ) < R 2 (Ai). Hence, there is spatial anti-monotonicity, directly implying 
spatial sub-additivity. All processes mentioned are stationary and have an isotropic and 
non-decreasing extremal coefficient function. For the Brown-Resnick model, we have, in 


the general case, 0(x, y) = 2$ 


7(*-y) 


, x, y e R 2 (see e.g. 


Davison et al. 


2012 ). 


3. First, we show that 


0 < / (2 — 0(x)) dx < oo yy 0 < a 2 < oo. 
Jr 2 


(29) 


Since u > 0 and, for all x e R 2 , 0(x) e [1,2], we have, using the mean value theorem, 
that 

0<exp(— gW),exp(— 2 )< 2 ~ 9(X) . 

V u J \ u J u 
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Hence, we immediately obtain that 


'E 2 


(2 — 0(x)) dx < oo =>- a 2 < oo. 


(30) 


Moreover, we have that 


0(x)\ f 2\ ( 2 A ( (2 — 0(x) 

exp l-— exp — = exp — exp - 

u J v uj V uj V V u 

>expf—0 2 -®< X > 


- 1 


u u 


which directly yields that 


a 2 < oo =>- / (2 — @(x)) dx < oo. 
J E 2 


(31) 


Moreover, since 


for all x G R 2 , 2 — @(x) > 0, j (2 — 0(x)) dx > 0 iff 2 — 0(x) > 0 for 

in 2 


all xfD, where PcR“ and satisfies \T>\ > 0. This directly gives 


'E 2 


(2 — 0(x)) dx > 0 a 2 > 0. 


(32) 


The combination of (30), (31) and (32) gives (29). Thus, (14) and (29) yield that 0 < 
a 2 < oo. Note that, although only the “=F’ implication of (29) is directly useful in this 


proof, we have also proven the implication in order to show that (14) is not more 
restrictive than the condition 0 < a 2 < oo. 

The remainder of the proof is partly inspired by the proof of Proposition 7.3 in 


Dedecker et al. (2007). Let A be a convex in A. For A > 0, we denote A\ = A A. 
For a set F C R 2 , let dF and Pp be the boundary and the perimeter of F, respectively. 


Let 


= A a , 


(33) 


for some a G (0,1). Let also d Vx A\ = {x G A\ : dist(x, dA\) < v\} and = 

7La\<% a 71a. Since A\ is compact, convex and has a positive Lebesgue measure, we easily 


obtain, using Theorem 2 (note that this result is related to Steiner’s formula) in Swanson 
(2011), that |<9„ A Ab\| = v\Pa x — 7rn 2 = \v\Pa — vrn 2 < Xv\Pa■ Hence, 


|<9„ a H a | ^ v x Pa 
A 2 |H| - A|H|’ 


which gives, using (33), that 


!“wpf = a < 34 > 

By definition of we obtain that = A 2 |H| — \d Vx A x \ 1 which gives, using (34), 

that 

|Af A) ' 


lim , 9I ,, 

A->oo \~\A\ 


= 1. 


(35) 
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have 


Using and the fact that, for all x,y G R 2 , 0(x,y) can be written @(x — y), 


we 


n 2 {\A) = 


A 4 |A^ 


IXAJXA . 


0(x — y)\ / 2 

exp l-— exp — 

u V u 


dx dy. (36) 


We introduce 


and 


k(x) = exp 


= 


©(x) 


A 2 |^| 


u 


IA x JAx 


U 


exp-, x G R 


k(x — y) dx dy, A > 0. 


Using (36) and (15), these notations lead to 


n 2 {XA) = -r^-rrTx and a 2 = / k(x) dx. 
X \A | J R 2 


(37) 


Note that, for all x G R 2 , 0(x) < 2, meaning that k is non-negative. We first show that 
lim^oo T\ = a 2 . We have, for all A > 0, that 


T\ — T\x + T 2) a + Ur a, 


(38) 


where 


Ti,x = 


To x = 


To x = 


-^ 2 |A| JA^ x) JyeA x -.\\x-y\\>v x 


k(x - y) dx dy, 


A |A| JA ( ” x) jyeA x .\\x~y\\<vx 


k(x — y) dx dy, 


A 2 |A| 


k(x — y) dx dy. 


'dv x A\ J A x 


Study of Ti x 

Using the facts that k is non-negative and that < A 2 |Al|, we obtain, for all 

A > 0, that 


\T h x\ < 


-^ 2 |A| JA ( ” x) J yeR 2 :|jx—y||>u A 


k{x - y) dx dy 


A 2 |A| . > A^ x ^ \./zeR 2 :||z||>'i;x 
< / k{ z) dz. 

■J z£R 2 :\\z\\>v\ 


k( z) dz I dx 


Thus, combining the fact that liniA^-oo v x — oo, the second part of (37) and the fact that 
a 2 is finite yields that 


lim Ti : x = 0. 

A—>-oo 


( 39 ) 


Study of T 2 ) 
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If x G ail( 4 y ^ 4a, we necessarily have that ||x — y|| > v\. Thus, we have, for 
all A > 0, that 


T 2l a = 


-^“|-4| JA^ x) •/yGR 2 :|jx—y|| <V\ 

i r 


^ 2 \A\ Ja ( : x) \J zSR 2 : 


k(-K - y) dx dy 
k( z) dz ) dx 


l|z||<«A 


1^4 


Ca)i 


A2|A| 


zeR 2 :|jz||<i,^ 


k( z) dz. 


Hence, we obtain, using (35), the fact that lim^oc V\ = oo and the second part of (37), 
that 


lim T 2 a = cr 2 . 

A—^oo 


(40) 


Study of T 3iA 

We have, using the fact that k is non-negative and the second part of (37), that, for 
all A > 0, 


|t 3 ,a| < 


^ 2 |"4| Jd v .A x J R 2 


k(x — y) dx dy = 


1 


^ 2 |"4| Jd v .A x \J R 2 


k( z) dz ) dx = 


Therefore, using (34), we obtain 


lim T 3ja = 0 


A—Kx) 


I^a^aI 2 

A 2 |H| 


(41) 


Finally, the combination of (38), (39), (40) and (41) gives that 

lim T x = cr 2 . 

A—>oo 


Thus, using the first part of (37), we obtain 


a 


A^|4|77. 2 (Av 4) ^ cr + o(l) i.e. ^(AA) 


+ ° 1 V 


Consequently, setting K\ = 0 and K 2 = tjt, the axiom of asymptotic spatial homogeneity 

1^41 

of order —2 is satisfied since cr 2 > 0. 

Let us now prove that the Smith process, the Brown-Resnick process with a semivar- 
iogram satisfying y(h) = ^||h|| a where rj > 0 and a G (0,2], and the tube process satisfy 
the required conditions. First, they are all simple and stationary. Thus, for all x, y G R 2 , 
@(x, y) only depends on x — y. Moreover, for all x G R 2 , 2 — @(x) > 0 on a set with 

positive Lebesgue measure. Hence, / (2 — 0(x)) dx > 0. 

Jr 2 

We now show that / (2 — 0(x)) dx < oo. In the case of the Smith process, recall that 
J R 2 


0(x) = 2$ , x G R 2 , (42) 

where ||x|| s = \/x'£ -1 x. Note that ||.||2 is the norm associated with the inner product 
induced by the matrix E _1 . Mill’s ratio gives us that the survival distribution function 
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$ of a standard Gaussian random variable satisfies <E>(/i) ~ 

h —^oo 


exp 


h?_ 
2 


(42), we obtain 


\Z2nh 


2 - 0 (x) = 2 


4> 


x s 


4 exp 


rv-/ 

II xll—>-oo 


v^ll 


x £ 


Hence, using 


(43) 


Since all norms in R 2 are equivalent, there exists B i > 0 such that, for all x G R 2 , 
||x|| s > Ri||x||. Thus, for all x G R 2 , 

4exp (-^) 4exp 


0 < 


v^H 


< 


B 2 ||x || 2 


x £ 


\Z2ttBi I 


(44) 


Using the fact that the right-hand term of (44) is integrable, the combination of (43) and 


(44) gives the integrability of 2 — 0. 


In the case of the Brown-Resnick process with a semivariogram satisfying 7 (h) = 7 ||h|| a 
where 7 > 0 and a G (0,2], the integrability is obtained in exactly the same way since 

0(x) = 2$ x e R 2 . 


In the case of the tube model, we have 0(x) = 2 for ||x|| > 2 R^. Therefore, the function 
2 — 0 has a compact support and is integrable. 


Finally, in all cases mentioned, we have / (2 — @(x)) dx < 00 . 

J R 2 


□ 


A.7 For Theorem [4] 

Proof. A random held {X(x )} xgR 2 is called associated if Co v(f(X I ),g(X I )) > 0 for 
any finite subset I C R 2 and for any bounded coordinatewise non-decreasing func¬ 
tions / : R card(/) R and g . R card(/) 1 —y R, where Xj = |A(x) : x G 1} and 
card stands for cardinality. A random held {X(x )} xgR 2 is called positively associated 
if Co v(f(Xi),g(Xj)) > 0 for all disjoint finite subsets I,J C R 2 and for any bounded 
coordinatewise non-decreasing functions / : R carc U / ) i_^ R and g ; R card(J) ^ R j£ jy j s 
associated, it is clear that X is positively associated. Proposition 5.29 in Resnick ( 1987j ) 
implies that max-stable processes are associated and therefore positively associated. 
Moreover, Z is stationary and measurable and thus continuous in probability. Since Z 
has standard Frechet margins, we obviously have that P(Z(0) = u) — 0 for all u > 0. 

For r > 0, we denote F +r = {x G R 2 : dist(x, F) < r}, where dist is the Euclidean 
distance. A Van Hove sequence in R 2 is a sequence (F n ) ne ^ of measurable subsets of 

I (dF n 


\+r I 


R 2 satisfying F n f R 2 , lim r 


FJ = 


00 , 


and lim,„ 


= 0 , for all r > 0 . 


For any convex A G A and any positive non-decreasing sequence (A n ) ng ]N such that 
lim^oo X n = 00 , the sequence (A n A) ng ]N is a Van Hove sequence. Indeed, for all n G IN, 
A n A C A n+ i A and (J^Li = R 2 , showing that the hrst condition is satished. The 
second one is trivial and we now show the last one. For all A > 0, since A\ is com¬ 
pact, convex and has a positive Lebesgue measure, we easily obtain, using Theorem 2 in 


Swanson (2011), that |(ch4A) +r | = 2 r Pa x = 2t\Pa- This directly gives 


lim 

A—Kx) 


\{9A, 


\ +r I 


A 2 |A| 


= 0 , 
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which immediately implies the last condition in the definition of a Van Hove sequence. 
Furthermore, it turns out that 


/r 2 


'r 2 


Cov (l{z(o)>u}> I{z(x)>u}) dx 


P(Z(0) > u, Z(x) > u) — 


'r 2 . 


exp 


0 (x) 


u 


exp 


u 


1 — exp — 
u 


dx = cr 2 . 


dx 


The combination of (16) and (30) implies that a 2 < oo. Moreover, > 0 since Z is 


measurable. Indeed, the case cr 2 = 0 corresponds to the situation where the random 
variables Z(x) are independent for almost all x G R 2 , which corresponds to processes 
which are not measurable. Therefore, we have that cr 2 G (0, oo) and, by applying Theorem 


7 in Spodarev (2014), we obtain, for all positive non-decreasing sequence (A n ) n ew such 
that limp 


A n = oo, that, for all convex A G A, 


\nV\A\ ( L]y(\ n A ) — m) —> A/"(0, a 2 ), i.e. X n ( L]\r(X n A ) — m) —> J\f 


a 2 

°'W 


(45) 


for n —>• oo. Now, let l : R —> R and B 2 G R. If, for any positive non-decreasing 
sequence (A n ) n£ ]N such that Hindoo X n = 00 , we have that lim^oo l{X n ) = B 2 , then 
lim^oo /(A) = B 2 . To prove this, let us assume that there exists a sequence (^ n )neiN such 
that lirrp^oc rj n = 00 but (/(?/„) ) ne in does not converge to B 2 . This means that there 
exists e > 0 such that, for all N G IN, there exists n > N such that |/(r/ n ) — B 2 \ > e. We 
introduce the sequence ( N n ) ne ^ defined by = inf{n G IN : 77 ^ > 0 and \l(r] n ) — B 2 \ > e} 
and, for n G IN, N n+1 = inf {ri G 1N\{1,..., N n } : \l{r] n ) — B 2 \ > e and rj n > rj Nn }• ft is 
clear that Ni and N n+ 1 exist since lim, woo r) n = 00 . By definition, the sequence (j 7 Ar n ) ne ]N 
is positive, non-decreasing, satisfies lini^oo r] Nn = 00 (as a subsequence of (rj n ) n£ i n) but 
does not converge to B 2 . This leads to a contradiction and shows the result. Hence, 
finally, we obtain 0 using ( ]45| ). 

Finally, the Smith process, the Brown-Resnick process with a semivariogram satisfying 
7 (h) = 7 ||h|| a where rj > 0 and a G (0,2], and the tube process are all stationary and 
measurable. Moreover, in all these cases, as shown in the proof of Theorem [3j ( Jl6| is 
satisfied. Thus, all the required conditions are satisfied. □ 


A.8 For Theorem [5] 

Proof. 1. The argument is the same as in the proof of Theorem [3j Bullet 1. 
2. By Theorem [4] and Proposition 0.1 in Resnick (1987), we have that 

a q a 


VaR a (A(L)v(AAl)) - m)) = X(VaR a {L N (XA)) - m) ->■ 


FT’ 


for A —> 00 . 


Thus, 

A(VaR a (LAr(AH)) — m) = + o(l) i.e. 77. 3 ,„(A A) = m + ~^-j= + o (\V 

A ^°° vAI A->-cxd Ay|H| \A ) 
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Hence, setting K\ = m and K 2 = —, the axiom of asymptotic spatial homogeneity of 

order —1 is satisfied since a > 0 and q a ^ 0 (because a ^ 0.5). 

Finally, regarding the examples of processes proposed, we know from Theorem [4] that 
the required conditions are satisfied. □ 
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